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SUMMARY 


1 . Background 

For  many  years,  the  Army  hat  used  progiam  factors  such  as  flying  hours, 
miles  travelled,  end  item  populations  and  rounds  fired  in  the  forecasting 
of  demand  for  repair  parts.  In  using  program  factors,  Che  assumption  la 
traditionally  made  that  demand  is  directly  proportional  to  the  level  of 
activity  - that  is,  an  increase  of  202  in  flying  hours,  for  example,  will 
cause  a 202  increase  in  parts  demand.  However,  a number  of  studies  have 
raised  doubt  about  this  assumption,  as  evidenced  in  reports  cited  in  the 
Selected  Bibliography  contained  in  thio  report.  Since  there  is  some 
conflict  in  the  evidence  of  these  studies,  this  work  was  undertaken  in  an 
attempt  to  resolve  the  issue. 

2 . Scope 

This  study  is  limited  to  developing  procedures  ,'nr  the  forecasting  of 
worldwide  recurring  demands  to'  Ati.y  managed  class  IX  secondary  items  (repair 
parts  and  spares)  including  stock  fund  and  PEMA  items.  Recurring  demand  ia 
that  portion  of  total  demand  which  is  due  ultimately  to  failures  of  parts 
in  service.  Specifically  excluded  from  recurring  demand  are  Initial  issues 
of  parts  when  new  units  are  activated  or  new  weapon  systems  are  fielded 
(provisioning),  issues  for  mobilization  reserves,  and  issues  to  maintenance 
facilities  for  rebuild  programs.  Also  excluded  are  those  secondary  items  that 
would  be  expected  to  receive  intensive  management  making  use  of  specialized 
Information  such  as  life-cycle  or  closed  loop  asset  data.  The  procedures  are 
to  be  applied  in  the  Commodity  Command  Standard  System  (CCSS)  inventory  man- 
agement function  implemented  at  the  Army's  National  Inventory  Control  Points. 

The  study  was  not  intentionally  limited  to  aircraft  items  but  the  only 
adequate  data  available  was  for  such  items-  However,  it  is  believed  that 
che  conclusions  may  be  extended  to  other  than  aircraft  items. 

3.  Methodology 

First,  a number  of  demand  models  were  considered  and  appropriate  fore- 
casting algorithms  were  selected.  These  algorithms  and  some  related  techniques 
were  then  screened  according  to  how  they  reacted  to  trends  in  aggregate  demands. 
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This  screening  was  done  to  weed  out  the  worst  algorithms  in  order  to  save 
resources  in  the  final  selection  phase.  Tl.e  ones  yielding  the  smallest  mean 
absolute  forecast  errors  were  then  compared  using  a simulation  of  the  in- 
dividual items  in  the  supply  system.  Tire  i Inal  selection  was  made  on  the 
basis  ut  smallest  aggregate  simulated  inventory  cost  for  constant  time-weighted 
requisitions  backordered. 

4.  Cone lusions 

The  Army's  presently  used  program  factor  is  better  than  a simple  moving 
average  of  demand  alone.  A forecast  method  based  on  a regression  equation  of 
demand  as  » function  of  flying  hours  is  better  than  Che  present  forecasting 
method . 

There  is  almost  no  improvement  to  be  gained  in  applying  different 
forecast  methods  to  different  segments  of  the  item  catalog  on  the  basis  of 
unit  price,  demand  rate,  dollar  value  of  demand,  or  dynamic  versus  nondynamic 
items. 

Finally,  there  is  no  improvement  in  program  factor  type  forecasts  due 
to  lagging  the  demand  behind  the  progiam  data. 

5 . Recommendations 

It  is  recommended  char  the  regression  algorithm  be  used  Instead  of  the 
currently  used  factor. 

Final  Implementation  should  be  delayed,  however,  pending  completion 
of  another  AMCIRC  study  on  Methodology  of  bemand  Forecasting,  which  is 
looking  at  methods  ocher  than  the  simple  moving  average  which,  when  combined 
with  the  program  factor  recommended  here,  night  give  better  forecasts. 
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CHAPTER 


PKr-  iOUS  WORK 

Previous  studies  of  the  relationship  between  demand  and  program  data 
by  the  Navy  and  Air  Force  have  met  with  1 ittle  success.  While  good 
correlation  has  been  found  for  gross  measures  of  demand  over  many  parts 
[Brown  1956,  DSA  196B]  the  relationship  ' etween  program  data  and  demand 
for  particular  parts  has  been  elusive  [ larger  1973,  Denicoff  1962, 

Gabriel  1972,  Haber  1967,  Mark! and  1970].  The  few  positive  findings 
do  not  relate  well  to  our  problem.  A study  by  RAND  tAstrachan  1961]  com- 
pared the  forecasting  accuracy  of  a simple  program  factor  to  that  of  three 
program  factors  with  built-in  safety  levels,  and  to  three  specialized 
methods  that  require  more  detailed  knowledge  of  part  history  than  i9  avail- 
able for  ordinary  repair  parts.  (These  service  life"  methods  are  applic- 
able to  expensive  components  like  aircraft  engines.)  Results  were  mixed, 
but  the  sample  program  factor  did  about  as  well  as  any  of  the  other  methods. 

The  authors  of  a thesis  presented  to  the  Air  Force  Institute  of 
Technology  [Goldfarh  1967]  stated  f’»at  program  data  had  "...little  value 
as  a demand  prediction  cool...";  however,  their  results  were  not  completely 
discouraging.  For  many  of  the  parts  studied  demand  correlated  well  with 
the  actual  flying  hours  of  tlu:  aircraft  » nJ  Items,  and  the  flying  hours 
themselves  could  bo  accurately  forecasted  from  the  planned  hours;  however, 
the  authors  were  cautious  about  drawing  conclusions. 

Previous  studies  have  used  forecast  error  as  the  criterion.  Thus  the 
best  forecast  method  v.  s considered  to  be  the  one  that  most  reduced  the  mea- 
sure of  deviation  betw<  en  forecasted  and  actual  demand  for  the  period.  How- 
ever, this  raises  the  question  of  which  . easurc  of  deviation  to  use.  Some  work 
has  used  mean  absolute  deviation  (MAD)  and  other  has  used  mean  squared  error. 
Results,  in  general,  may  be  different.  Another  consideration  is  that  some 
items  have  a larger  impact  on  supply  performance  than  others  - and  that  some 
items  cost  much  more  than  others.  For  example,  a study  with  Navy  data  [Denicoff 
1962]  concluded  "...only  59  parts,  or  three  percent  of  the  parts  examined, 
were  correlated  with  flying  hours,''  The  author  did  not  know  whether  this 

Best  Available  Copy 
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group  of  parts  represented  3 percent  or  30  percent  of  the  total  parts 
budget.  It  was  for  this  reason  that  a cost  versus  performance  simulation 
approach  was  prescribed  by  the  Office  of  the  Secretary  of  Defense  for 
Installations  and  Logistics  as  the  method  to  be  used  in  validating  new 
forecasting  techniques It  has  already  been  used  by  Inventory  Research 
Office  to  develop  the  details  of  the  inventory  policy.  See  "Evaluation 
of  Several  V5L/E0Q  Models"  by  Robert  L.  Deemer  and  W.  Karl  Kruse,  May 
1974,  AD  781948,  "Frequency  of  Requirements  Determination"  by  Robert 
L.  Deemer,  October  1974,  AD  A003227,  and  "Estimation  of  Demand  Variability 
Parameters"  by  Alan  J.  Kaplan,  May  1974,  AD  781942. 
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2.1  Source  and  Organization 

The  conclusions  in  this  report  a:o  ! ased  on  studies  made  using 
chronological  demand  data  from  AVSCOM.  Hying  hours  covering  the  same 
period  as  the  demand  data  were  obtained  from  PCSLOG.  All  data  has  been 
summarized  by  quarter.  For  each  quarter  we  have  the  worldwide  totals  of 
the  number  of  requisitions,  the  quantity  detaar.ded  and  the  haurs. 

The  Hying  hour  totals  are  broken  out  by  aircraft  type /model /series  (TMS) . 
The  data  spans  the  28  quarters  from  Janu.arv  1967  thru  December  1973. 

2.2  Weapon  Systems 

Appendix  A,  Section  V or  AF  710- i relates  end  items  into  classes 
called  weapon  systems.  A two  letter  cod*  is  assigned  to  each  system. 

The  first  letter  is  an  A for  fixed-wing  ilrplanes,  and  B for  helicopters. 
The  table  below  contains  the  TMS  codes  gieuped  into  weapon  systems. 
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A chart  on  the  next  page  gives  an  idea  o the  appearance  of  the  various 
aircraft.  The  heading  "DESIGNATION"  corresponds  to  TMS. 

The  graph  on  the  following  page  and  those  in  Appendix  2 6how  the 
flying  hours  and  demand  quantity  by  quarter  lor  the  first  five  years 
for  those  weapon  systems  with  a reasonably  large  group  of  repair  parts. 

The  demand  and  hours  have  been  scaled  (as  described  in  the  chapter  on 
Model  Valuation  and  Selection)  so  that  each  sums  up  to  one  (unit  demanded 
or  hour  fl >wn)  for  the  20-quarter  period  for  each  part.  The  graphs 
show  the  a/erage  of  the  scaled  values  by  quarter  for  a large  sample 
of  the  parts  in  each  system  (the  preliminary  data  base  described  below). 

It  can  be  seen  from  the  graphs  that  the  demand  follows  the  flying 
hour  pattern  fairly  well  on  some  of  the  systems.  The  exceptions  are  mostly 
in  the  form  of  peaks  in  demand. 

These  observations  suggest  that  demand  is  well  correlated  to  flying 
hours,  and  that  a program  factor  would  be  expected  to  yield  good  demand 
forecasts  by  the  traditional  forecast  error  criterion,  at  least  for  this 
data. 
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HELICOPTER  SERIES 


OESIGN’AliOM 
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AGCRECATE  NORMALIZED  DEMAND  QUANTITY  AND 
FLYINC  HOURS  FOR  89  U-21  PARTS 


Legend : 

F » Flying  Lours 
D “ Demand  Quantity 
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Note: 


See  Appendix  2 for  additional  aircraft. 

10 


2 . 3 Input  Data  Editing 


If  an  item  was  used  on  more  than  on1  aircraft  TMS  in  a system,  the 
quantity  used  per  aircraft  and  the  venrout  rate  were  assumed  to  be  similar 

on  each  aircraft.  Consequently,  for  these  items,  the  flying  hours  for  all 

applicable  aircraft  in  the  system  were  simply  added  together.  Parts  applying 
to  TMS  models  of  different  weapon  systems  were  not  included  in  this  study. 

Flying  hour  data  were  available  only  for  US  Army  and  Vietnam  Air 

Force  programs.  Because  of  this  limitation,  processing  was  begun  on  a 

(voluminous)  demand-transaction  file  with  customers  identified  in  order 

to  remove  demands  from  other  military  services  and  MAP  countries.  Some  of 

the  study  results  are  based  on  a preliminary  data  base  developed  in  this 

way  from  part  of  the  detailed  demand-transaction  file.  The  eleven  weapon 

systems  each  with  acceptable  data  for  at  least  40  parts  are  marked  with  a 

dash  in  the  weapon  system  table  above.  uere  are  about  2000  items  in  this 

preliminary  data  base.  It  covers  the  5-year  period  from  January  1967 

thru  December  1971,  The  data  were  aggregated  by  weapon  system  as  shown 

in  the  graphs.  In  developing  the  preliminary  data  base  it  was  discovered 

/ 

that  the  non-Army  demands  were  a small  fraction  of  the  total  number.  So 
the  main  part  of  the  work  was  done  with  n larger  data  base  which  had 
already  been  aggregated  to  include  all  customers.  (Only  the  customers 
were  aggregated  - not  the  items.)  The  advantages  to  be  gained  in  using 
the  larger  data  base  were  great  enough  to  warrant  the  small  risk  that  for 
a few  individual  items  the  non-Army  demands  might  actually  constitute  the 
major  part  of  their  demand.  The.  following  table  gives  the  number  of 
parts  (total  11631)  by  weapon  system  for  the  expanded  data  base. 


DISTRIBUTION  OF  PARTS  BY  WEAPON  SYSTEM  IN  AVSCOM  7-YEAR  DATA  BASE 
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The  scope.  of  tills  work  war,  .limited  to  recurring  demand.  When  the 
demand  waa  accumulated  Into  quarters  the  following  classes  of  requisition 
were  not  included: 

initial  issue, 
mobilization  and 
rebuild 

If  the  remaining  demand  was  so  small  that  there  was  not  even  one 
calendar  year  with  at  least  four  requisitions,  the  item  was  dropped.  Such 
items  are  called  trivial-demand  items  belov:. 

The  rest  of  this  chapter  applies  only  to  the  large  data  base  used  for 
the.  main  part  of  the  study.  If  there  were  no  flying  hours  in  the  first  two 
quarters  (January  to  June  1967)  the  item  was  considered  to  bo  new  to  the 
system  and  was  dropped  (because  it  would  likely  be  managed  with  maintenance 
factors  in  transition  from  provisioning.)  likewise  an  item  for  which 
there  were  no  flying  hours  in  the  final  two  quarters  (July  to  December  1973) 
was  not  included  (because  the  item  manager  would  have  tempered  any  dev, land 
forecast  with  the  knowledge  that  the  airctaft  was  being  phased  out).  Finally, 
items  with  over  a million  dollars  cf  demand  per  year  - see  below  for  how 
estimated  - were  dropped  because  there  arc-  special  forecasting  methods  for 
such  parts.  Each  itetji  was  classified  as  low  or  medium  dollar  value  (LDV) 
or  high  dollar  value  (HDV)  according  to  whether  the  demand  rate  averaged 
over  the  full  28  quarters  was  less  than  or  greater  than  or  equal  to 
$50,000;  and  the  requisition  rate  was  less  than  or  greater  than  or  equal 
to  100  per  year. 

To  summarize,  items  were  excluded  for  the  following  reasons: 

38%  managed  by  DISC  or  other  JUT  or  not  stocked 
28%  trivial  demand 

19%  no  way  to  associate  with  flying  hours 
2%  applied  to  more  than  one  system 
1%  new,  phase-out  or  super  high,  dollar  items 

This  information  is  given  in  more  detail  in  the  "pipe"  chart.  Near 

the  bottom  of  the  chart  is  the  legend  "lOUT  GOOD  ITEMS  LEFT  T.N  THE  PIPE 

AT  THIS  POINT".  The  additional  items  removed  below  that  point  were  dropped 

for  two  reasons.  First,  half  of  the  I.DV  non-dynamic  items  (defined  later)  - 

including  18  Pi"./,  and  5137  stock  fund  - were  dropped  to  reduce  computer 
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DATA  EDITING 


Pipe  Contains 
88944  Items 
in  AVSCOM 
Catalog 
Dec  1971 


33995  Items  Not  Managed  by  AVSCOM 
or  Not  Repair  Parts 


24717  Items  With  Trivial  Demand  in  1967-1973 


16727  Items  For  Which  No  Program  Data  Available 


1835  Items  Applied  to  More  Than  One  Weapon  System 


- 


..  v 


1122  Items  New  in  1967  or  Phase  Out  By  1973 


39  Items  With  Annua  l._P.emand  of  More  Than  $1  Million 

10509  Good  Items  Left  in  the  Pipe  at  This  Point 
JT  5186.  ADV  Non-Dynamic  Items  Not  Used  (Including  29  PEMA) 

Pipe  to  Simulator  Now  Contains  the  Remaining  5323  Items 


expense  in  the  simulation  runs.  To  compensate,  we  double  weight  the  LDV  non- 
dynamic results  in  the  final  tabulations.  Then  the  remaining  11  PEMA  items 
were  also  dropped.  The  reason  this  was  done  was  because  the  PEMA  items  were 
such  a tiny  fraction  of  the  I,1)V  non-dynamics  (?.9  out  of  10350)  that  they 
were  thought  to  be  atypical.  Usually  PEMA  items  are  expensive  and  repar- 
able items;  the  LDV  non-dynamic  class  was  intended  to  contain  relatively 
cheap  non-reparablcs. 


I 

I 

2.4  Final  Characteristics  of  Data  Base 

The  data  base  was  separated  ii.u-  l.l>v  and  UDV  portions  as  described  1 

above.  Kach  of  the  portions  was  divided  further  into  dynamic  and  j 

non-dynamic  portions.  The  criterion  for  this  partition  was  the  Item's  j 

Federal  Stock  Class  (FSC) . The  table  below  Rives  the  fraction  of  items 
in  thu  most  important  aircraft  Federal  Stock  Classes.  The  upper  part 
of  the  table  contains  the  classes  included  in  the  definition  of  c/namic 
items.  The  important  dynamic  components  were  considered  to  be  those  that 
experience  very  high  rotation  rates;  i.e.,  rotor  blades,  transmissions,  i 

turbine  * rgines  and  other  engine  components.  This  encompassed  about  ten 
percent  oi  the  items  (FSCs  selected  for  us  by  AVSCOM  to  contain  dynamic  j 

items  are  1615,  2840  and  2810).  ; 


DISTRIBUTION 

OF  ITEM  FSC  CLASS 

PERCENT 

FED.  STK  CLASS 

DESCRIPTION  (AR  708-15) 

9 

1615 

Helicopter  Rotor  Blades  and  Drive 
Me  haniams  (Transmissions) 

i 

2810/2840 

Engines 

10 

TOTAL 

Dynamic  Items 

! 

i 


39 

1560 

Airframe  Structure 

14 

1680 

Misc.  Accessories 

4 

1650 

Hydraulic,  Vacuum  and  De-icing 

3 

2995 

Engine  Accessories 

2 

3120 

Bearings,  Plain,  Unmounted 

2 

5340 

Misc.  Hardware 

2 

1620 

Landing  Gear 

2 

1660 

Air  Conditioning,  Heating  and  Pressurizing 

1 

5330 

Packing  and  Gasket  Materials 

21 

other 

90 

TOTAL 

Non-E  namic  Items 

15 


Notice  that  the  nondynamic  items  include  bearings  and  engine  accessories 
some  of  which  are  likely  to  have  dynamic  characteristics.  Thus  this  class 
may  be  contaminated.  It  was  felt  ’.'.at  i«.  was  more  important  to  have 
an  uncontaminated  group  of  dynamic  items  than  nondynamic  in  order  to  be  able 
to  draw  certain  conclusions  about  forecasting  methods. 


The  final  data  base  characteristics  are  shown  in  the  following  table.  The 
major  classes  of  items,  with  subtotals  ate  given  in  the  left  hand  column. 
The  next  column  gives  the  total  number  oi  items  in  each  class.  The  reason 
the  total  is  11631  here  (and  in  the  weapon  system  chart)  Instead  of  10509 
as  in  the  "pipe"  chart  is  because  this  table  uai  aide  before  the  1122  NEW 
AND  PHASE-OUT  items  lud  been  removed.  The  las.  two  columns  give  the  totals 
broken  out  by  fund  (PEMA  or  ASF)  and  segment  (non-reparable,  reparable 
or  insurance  item) . 

DATA  HASP.  CHAKACT'  RISTICS 
DISTRIBUTION  BY  ITEM  < HARACTERISTICS 
FINAL  ITEM  SET  IN  IMULATION 


ITEM  CLASS 

TOTAL  ITEM  COUNT 

PEMA 

ASF 

NON-REP 

REP 

INS 

LDV  Non-Dynamic 

10350 

29 

10321 

! 9877 

448 

25 

LDV  Dynamic 

1008 

24 

984 

957 

51 

LDV  TOTAL 

11358 

53 

11305 

10834 

499 

25 

HDV  Nor -Dynamic 

174 

30 

144 

66 

108 

HDV  Dynamic 

99 

52 

47 

26 

73 

HDV  TOTAL 

273 

82 

191 

92 

3.81 

TOTAL  Non-Dynamic 

10524 

59 

10465 

9943 

556 

25 

TOTAL  Dynamic 

1107 

76 

1031 

983 

: :>4 

TOTAL  Cood  Items 

11631 

135 

11496 

10926 

680 

25 
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CHAPTER  III 


MODEL  DEVELOPMENT 

We  consider  three  simple  alternative  demand  models: 

1.  Average  demand  is  constant. 

2.  Average  demand  is  proportional  to  total  flying  hours. 

3.  Average  demand  has  a component  described  by  Model  1 and  a 

component  described  by  Model  2. 

In  none  of  those  three  models  is  the  demand  to  be  any  different  from  failures 
ce  parts  at:  the  user  level.  However,  in  the  fourth  model  below,  we  implicitly 
allow  for  the  delay  due  to  stockage  at  supply  levels  between  the  original 
failure  and  the  final  demand  at  the  NICP  level. 

•'*.  Failures  are  proportional  to  total  operating  hours;  a demand 
is  delayed  from  the  time  of  the  failure-  by  an  unknown  amount. 

In  all  four  models  the  word  "constant  *’  means  a value  that  may  change 
gradually  but  does  not  depend  on  demand  or  program.  The  forecasting 
algorithms  to  be  tested  are  based  on  these  models.  They  are  grouped  under 
the  appropriate  model  following  the  definition  of  the  notation. 

Notation 

Subscripts  represent  time  as  follow;-:: 

Small  letters  stand  for  quarters,  largo,  letters  stand  for  several 
contiguous  quarters.  Time  periods  are  always  relative  to  the  present;  the 
forecast  is  being  computed  as  of  the  end  of  quarter  0.  There  are  B quarters 
in  the  base  period  and  P quarters  in  the  forecast  period. 

Base  period: 

where  D is  the  actual 

q 

demand  in  the  qth  quarter 

where  F^  is  the  program  data 

value  (flying  hours)  in  the 
qth  quarter 


o 

n-  i d 

q—  (B-l)  q 


o 

FB  * F F 


q—  (B-l)  q 
✓ 
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Forecast  period: 


AUP 


P 


v 

q* 


1 


r 

q 


This  1e  the  final  forecast  for  algo  Ithm  A.  I)  is  tlu-  estim-ite  given 

q 

by  algorithm  A far  the  demand  in  the  qth  quarter. 

In  the  following  descriptions  of  tl\  - ulgoritluas  any  unspeciliec  suns 
are  intended  to  be  over  the  base  period. 

Model  1 

Demand  is  constant. 

yjVAVo  Algorithm  (moving  -tv<a  ige) 

"p  ■ <t>  "b 


Comment : 

The  "constant"  demand  rate  is  stimated  as  the  mean  of  rcient 
observations. 
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TKF.XPF  Algorithm  (triple  exponential  sn  >thlng) 


I) 

P 


P 

II 


q-1 


1) 

q 


where 


’ flt  + btq  + T Sq2 


ac  - Do  + U-Y)3  (D^  - Po) 


fct  * Vi + ct-r*^(2-Y)(,)t-i 


V 


CC  ‘ Ct-1  ‘ ^ <UUl  ' V 


Dc  “ at  + bt  + 2 Ct 


and  y is  the  fixed  smoothing  constant. 

Comment : 

The  t index  represents  the  absolute  quarter;  the  q index  is 
for  the  quarter  relative  to  the  present.  Thus  Pq  is  the  actual  demand  in 
the  quarter  just  ended. 

The  a,  b,  c initial  values  are  developed  from  a least  square 

linear  regression  on  the  3 earliest  available  quarters  of  demand.  The 

initial  value  of  P is  set  equal  to  this  initial  a-value  unless  an  estimate  is 
o 

available. 

The  algorithm  u is  tried  with  the  following  values  of  y to  give 
four  separate  estimates  of  demand: 


Y • .05,  .1,  .15.  .19 

In  this  algorithm  the  gradual  changes  in  th:  "constant"  demand  rate  are 
given  explicit  recognition  ir  the  trend  istimators  b and  c. 

Reference:  Robert  Goodell  Brown,  Smoothing  Forecasting  and  Prediction  of 

Discrete  Time  Series,  Prentice-llall  1963,  Chapter  9. 
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Model  2 


Demand  is  proportional  to  total  flying  hours. 
SEV94  Algorithm  (Issue  Interval) 


n 


p 


£>n 

s 


a 


C onnnent : 

. this,  torni  It  is  viewed  as  a moving  average  correct  h\  a 
program  (actor.  It  tray  also  be  written: 


BP  • (£  F,, 


In  tills  equivalent  form  it  is  seen  to  he  a relationship  between  uii.and  anti 
program  data  like  the  Following  algorithm. 

REGORG  Algorithm  (Slope  linear  regression) 


where 


in  f 

b-— 4A 


I F 


Comment : 

Notice  that  both  the  l.EGORG  and  SEV94  algorithms  arc  spi  .iai 
cases  of  the  form: 

Iw(q)D 

D ■ ■ - p 

p E«(q)Fq  p 


where  tnc  weight  function  w .s  constant  in  the  SEV94,  but  is  variable 
in  the  REGORG  and  equal  to  F . Thus,  it.  tht?  SEV94  algorithm  eaci  quarter 
is  given  weight,  in  the  REGORG  algorithr  each  flying  hour  is  given  equal 
weight . 
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« 


f 

i 

l 


% 


SBAS ) .:  1 Algorithm  (sliding  base  period) 


l) 


1' 


!> 

f i*  ’ 


i)' 

q 


q-i 

r 

L 

iia i 

q-l 


I F 
i-q-H 


i 


F 

q 


In'  if  i ••  o 
D}  i f i •_  o 

Comment : 

This  algorithm  is  similar  to  the  SKV94  algorithm..  The  diflirence 
is  that  each  successive  quartet  is  the  forecast  period  is  forecast  using 
a base  period  which  ends  with  the  previous  quarter.  For  all  except  the 
first  quarter  this  means  that  forecasted  demand  appears  in  this  base  period. 
This  results  in  more  weight  being  given  to  recent  quarters  than  to  earlier 
quarters  of  the  real  base  period  in  any  forecast  for  more  than  one  quarter 
in  the  future. 


2J 


im-imiihi*.  II*  MfUa I > 


Model  3 


Demand  is  partly  constant  and  partly  proportional  to  operating  hours 
PARTEF  Algorithm  (partial  effect) 

Dp  * ^Y)(SEV9ADp)  * (1  ^ ^AVGMOV° p^ 

where  y Is  a constant. 

Comment : 

This  is  a weighted  mean  of  two  simple  estimates  described  above 
Y is  intended  to  be  fixed,  but  possibly  re-evaluated  from  time 
to  time,  and  will  be  the  same  for  all  items. 

The  algorithm  was  tested  with  the  following  values  to  give  four 
separate  estimates  of  demand: 

y « .25,  .5,  .75,  ,95  , : _ ~“ 
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■ 


i 


m 


THkKSll  i tj_nii  (program  lactor  only M_  n discriminant  function  reaches 

w«e-  a threshold  value) 


Dp  " ^R)(RKGORcV  + (1'L;(AVGMOVDP) 

where 

B is  0 or  1 as  deCerrnlned  by  the  following  test: 
b 

Let  t * 


s iu  an  estimate  of  the  standard  deviation  of  the  error  In  the 
fitted  line: 


l> 

q 


b F 

q 


T 

BD"  - 

q 


id2f 

. <LJL 

IF2 


df 


and  df  » Base  Period  - 1. 

if  t > tdf 


Then 


B‘t 


otherwise 


(2.35  if  df  - 3 

1.89  if  df  ■ 7 

1.80  if  df  = 11 


Comment : 

This  algorithm  Is  a heuristic  that  determines  whether  the  flying 
hours  and  demand  are  closely-cnough  related  in  the  base  period  to  warrant 
using  the  REGORC  forecast.  If  t is  too  :mall  the  algorithm  falls  back 
on  the  noving  average. 
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**  > r-  ' 


PARTRC  Algorithm  (PARTEF  with  variable  parameter) 


DP  “ b ^SKV^P*  + (1"b)  (AVCMOvV 


where  I)  Jh  vat iuble 


Error 


V - V 

max  (l,I)pt 


where  D_,  “ I D 

P q-(P’-l)  q 


and  Dp,  ■ bOgr)  + p 

B 


and  B'  - l*'  - -7 


Coimnent : 


The  parameter  h is  determined  as  the  value  in  range  [0,  .05, .1, . • .1.0] 
which  minimizes  the  error  in  forecasting  the  second  half  of  the  base  period 
using  the  first  half  as  a base.  If  the  error  > .4  b is  set  to  0. 

KEG I NT  Algorithm  (Slope,  Intercept  linear  regression) 


Dp  - a + b Fp 

p 

where  a ■ (ij)  " b *’g) 


Fr 

l (F  - ~ -)  (U  - "T”) 
q B Q B 


Comment : 


Note  that  the  PARTEF  algorithm  reduces  to  the  RECINT  if  vp  make 


the  following  subs -itut ion: 


y = b instead  ci  using  a fixed  value. 


POWKRF  Algorithm  (fixed  power) 


y 

AVCMOV  UP 


0 < y < 1 


where  y is  fixed. 

Comment : 

In  POWERF,  y is  used  in  the  sam  way  that  it  is  in  the  PARTEF 
algorithm.  This  algorithm  has  a compromise  character  like  the 
rARTF.F , but  its  form  is  non~ilnear.  The  parameter  y determines  the  fraction  of 
the  rate  of  change  of  flying  hours  by  which  the  demand  forecast  will  change. 
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I 


POWREC  Algorithm  (variable  power) 


D 


P 


when: 


D 

q 


b 

P (AVGMOVUp) 


and 


b 


■~y  constrained  to  - 10^  <_  b <_  1 

zr 


b 

where  Y - log 

B 

n 


Y 


D 


Comment : 

This  algorithm  was  developed  by  linear  regression  on  the  logarithm 
of  the  power  form.  Note,  in  the  formula  for  b,  if  either  X or  Y is  undefined 
for  a quarter,  that  quarter  is  excluded  from  the  sums  in  both  numerator  and 
denominator. 
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Model  4 


Failures  are  proportional  to  total  operating  hours;  a demand  is 
delayed  from  the  time  of  the  failure  to  uliov  for  the  effect  of  stockage 
at  intermediate  supply  points. 

To  test  this  model,  each  of  the  algorithms  that  made  use  of  program 
data  were  tested  with 


F replaced  by  F „ 

q q-K 


where  the  three  values  1,  2 and  3 quarters  were  tried  for  the  lag  K. 
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EXAMPLE  OF  APPLICATION  OF  SOME  QI  THE  FORECAST  ALGORITHMS 


HYPOTHETICAL  DEMAND  AND  FLYING  HOUR  DATA 


QUARTER 

DEMAND 

FLYING  HOURS 

Jan  67 

12 

230 

Sep  67 

0 

245 

Dec  67 

0 

271 

Mar  68 

3 

270 

jtm  68 

1 

250 

Sep  68 

5 

220 

Dec  68 

0 

200 

Mar  69 

8 

190 

END  OF 

2-YR  BASE  PERIOD 

Jun  69 

185 

Sep  69 

187 

Dec  69 

180 

Mar  70 

END  OF 

1-YR  FORECAST 

PERIOD 

160 

FORECASTED  DEMAND 

ALGORITHM 

FORECAST 

ALGORITHM 

FORECAST 

MOVAVG  & PARTRG 

14.50 

REGINT 

27.34 

SEV94 

11.01 

S BAS EM 

9.40 

PARTEF.5 

12.  -.3 

POWREG 

51.10 

THRESH  & KEGORG 

10.27 

TREXP.05 

1324.73 

P0WERF.5 

12.63 
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CHAPTER  IV 


MODEL  VALIDATION  AND  SELECTIOH 

4.1  Screening 

The  method  by  which  the  final  selection  of  the  best  algorithm  was  made 
is  described  in  the  next  section.  In  this  section  is  detailed  how  the 
algorithms  described  in  the  preceding  chapter  were  first  screened  to  reduce 
them  to  a small  enough  number  so  as  to  be  economically  compared  by  simulation. 

It  was  originally  supposed  that  most  of  the  algorithms  would  do  poorly 
compared  to  a few  good  ones.  However,  as  shown  in  the  chapter  on  Results, 
most  of  the  algorithms  performed  tbout  the  same  in  screening. 

The  screening  criterion  wes  sufficiently  different  from  that  of  the 
simulation  to  give  additional  confidence  in  an  algorithm  which  did  well  in 
both.  Since  many  algorithms  performed  similarly  in  screening,  the  selection 
of  algorithms  for  final  comparslon  was  somewhat  subjective.  The  moving 
average  algorithm  (MOVAVC)  was  Included  as  a control  in  order  to  determine 
whether  the  use  of  a program  factor  had  any  value  at  all.  (Since  aggregate 
demand  followed  the  program  data  so  closely,  as  can  be  seen  in  the  graphs 
in  the  Data  chapter,  it  seemed  likely  th.it  individual  item  demand  would 
also. ) 

In  addition,  the  Army's  presently  used  forecasting  algorithm  (SEV94) 
had  to  be  included  in  the  final  testing  so  that  its  relative  value,  as 
compared  to  other  algorithms,  could  be  assessed. 

It  was  decided  to  include  only  two  more  algorithms  so  that  the  avail- 
able computer  funds  could  be  spent  in  more  Intensive  testing,  l.e.,  more 
Items  In  the  data  base. 

The  screening  process  will  now  be  described  in  detail.  The  data 
was  first  aggregated  by  weapon  system.  Equal  weight  was  given  to  each 
item  independently  of  its  cost.  This  was  done  by  normalizing  the  demand 
so  that  the  relative  quarter-to-quarter  variations  counted  equally  for 
different  items.  The  normalized  velue  fur  quarter  q is: 
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for  each  Item. 


D 

. -a._ 

20 

>'•  D 

q-1  q 


Then,  this  quarterly  value  was  averaged  over  all  the  Items  In  the 
weapon  system.  The  example  should  make  this  normalizing  process  clear. 


NORMALIZATION  EXAMPLE 

(FOR  CLARITY  6-MONTH  PERIODS  ARE  L1IOWN  INSTEAD  OF  QUARTERS) 


YEAR 

1967 

1968 

1969 

1970 

1971 

TOTAL 

ITEM 

0 

0 

1 

1 

5 

4 

O 

5 

0 

20 

1 

0 

0 

.05 

.05 

.25 

• *- 

t 

• A. 

.1 

.25 

0 

1 

ITEM 

0 

0 

6 

6 

30 

24 

12 

12 

30 

0 

120 

2 

0 

0 

.05 

.05 

.25 

• •- 

.1 

.1 

.25 

0 

1 

ITEM 

0 

4 

3 

0 

0 

1 

0 

0 

2 

0 

10 

5 

0 

.4 

.3 

0 

0 

.1 

0 

0 

■y 

* — 

0 

1 

NORMALIZED 

AGGREGATE 

0 

.133 

.133 

.033 

.167 

.167 

.067 

.067 

.233 

G 

FOR  EACH  ITEM  THE  UFPER  LINE  IS  THF.  ACTUAL  DEMAND, 
THE  LOWER  IS  THE  NORMALIZED  VALUE 


In  the  example,  items  1 and  2 have  exactly  the  same  quarter-to-quarter 
variations  even  though  their  demand  totals  are  quite  different.  The  third 
item  shows  much  more  variation  chan  the  first  two.  Notice  that  the  final 
aggregate  tends  to  be  smoother  than  the  individual  items.  Because  of  this 
an  algorithm  that  does  well  in  the  screening  might  do  poorly  in  the  final 
selection  when  individual  items  are  used.  This  should  be  especially  true 
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of  an  algorithm  like  triple  exponential  smoothing,  which  responds 
radically  to  sudden  changes. 

The  forecast-error  criterion  Is  used  merely  as  a screening  method. 
Final  selection  Is  based  on  economic  considerations  to  be  described  below. 
The  measure  of  forecast  error  employed  in  the  screening  phase  Is  the  mean 
absolute  value  (MAD).  The  MAD  was  used  rather  than  a square*,  measure  so 
as  not  to  overemphasize  large  errors. 

The  testing  procedure  was  designed  to  yield  the  greatest  number  of 
forecasts  while  avoiding  overlapped  forecast  periods.  An  example  will 
help  to  explain  this  procedure.  Assume  we  are  testing  with  a two- 
year  base  period  and  a one-year  forecast  period.  The  base  period  f^r 
the  first  forecast  will  cover  the  first  two  years,  and  the  first  forecast 
period  will  Immediately  follow  the  end  of  this  base  period.  After  the 
forecast  is  compared  to  the  actual  demand  in  the  forecast  period,  and 
the  resulting  error  saved,  the  process  is  repeated  for  the  next  forecast 
period.  To  do  this  the  base  and  forecast  periods  are  shifted  forward  by 
the  length  of  the  foreca-t  period  (one  year  in  the  example). 
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. ... 

FORECAST 

** 

SECOND 
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1967 

1968 

1969 

1970 

r 1 
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As  can  be  seen  there  Are  three  forecasts,  and  thus  three  values  of  MAD, 
for  the  case  of  a two-year  base  period  and  a one-year  forecast  period. 

For  screening,  three  different  forecast  period  lengths  (.5,  1 and  1.5  years) 
were  employed  and  the  results  averaged.  The  table  gives  the  number  of  MADs 
for  each  forecast  period  length  for  each  of  two  base  period  lengths  studied. 


TOTAL  FORECAST  MADs  FOR  EACH  WEAPON  SYSTEM 


BASE 

FORECAST  PERIOD 

TOTAL  FOR 

PERIOD 

.5  YRS 

1 YEAR 

1.5  YRS 

BASE  PEFvIOD 

2 YRS 

6 

3 

2 

11 

3 YRS 

4 

2 
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In  comparing  the  algorithms  for  a 2-year  base  all  eleven  forecast 
raaulta  could  be  used.  However,  to  compare  base  periods,  only  7 could 
be  used  because,  to  be  strictly  comparable,  the  forecasts  must  cover  the 
n m rime  period;  and  the  year  1969  is  excluded  as  a forecast  for  the 
3-year  bar-'  period.  A similar  remark  applies  to  testing  the  etfects  of 
lag;  in  this  case  nine  forecasts  were  available  since  the  maximum  lag 
was  2 quarters. 

A. 2 Simulator  Overview 

The  final  selection  among  algorithms  was  made  by  observing  their 
operation  in  the  Army's  NICP  context  using  the  DoDI  4140.39  simulator. 

This  simulator  is  described  in  the  report  "ALPHA  4140.39  Simulator"  by 
Martin  Cohen  (May  1973)  AD  762348.  Some  changes  were  made  to  the  simulator 
since  the  report  was  published;  the  major  changes  are  listed  in 
Appendix  1. 

The  general  organization  of  the  simulator  is  given  in  the  chart,  which 
shows  the  flow  of  information  among  the  processes. 
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The  demand  history  for  a given  item  is  fed  both  to  the  forecast  algorithm 
process,  which  computes  the  demand  forecist  for  each  period;  and  to  the 
asset  monitor,  which  reduces  on-hand  stock  - thus  possibly  triggering 
an  order  for  replenishment  - and  computes  the  final  performance  estimates. 
The  demand  forecasts  are  fed  to  the  inventory  model,  which  computes  the 
reorder  point  needeu  for  asset  monitoring  end  the  order  quant  Ity  needed 
when  orders  are  placed.  The  final  performance  estimates  are  averaged 
over  the  simulated  time  for  each  Item,  and  then  aggregated  across  all 
items.  This  is  done  separately  for  each  forecasting  algorithm.  The 
resulting  values  are  compared  to  determine  the  best  one. 

The  results  of  the  simulation  runs  appear  in  the  next  chapter  in 
the  form  of  graphs  of  performance  versus  cost.  A typical  graph  looks 
like  this. 


The  curves  are  traced  thru  several  points  for  each  algorithm;  each  point 
results  from  a run  of  the  simulator  with  another  value  of  the  lambda  para- 
meter. (See  the  Simulator  report  or  the  report  by  Deemer  and  Kruse  cited 
In  Appendix  1.)  In  the  above  hypothetical  graph,  Algorithm  A is  seen  to 
be  superior  to  Algorithm  B because  it  incurs  a lover  cost  for  the  same 
performance. 
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CHAPTER  V 


RESULTS 


5.1  Algorithm  Selection 

Th«  results  of  screening  arc  shown  on  the  bar  graph.  Recall  from 
section  4.1  that  the  criterion  for  comparing  algorithms  is  mean  absolute 
deviation  of  the  forecast  from  the  actual  demand  (MAD) . This  MAD  is  the 
mean  of  121  samples  (two-year  base  period  with  forecasts  for  6 months,  one 
year  and  18  months  - 11  forecasts  for  each  ol  11  aircraft  systems).  The  best 
algorithm  has  a MAD  of  .043;  since  the  expected  value  of  the  forecast  is 
.16,  thi9  is  a relative  error  (PCER)  of  26%. 

Notice  that  the  difference  between  the  best  and  worst  algorithms 
that  utilize  program  data  is  about  the  same  as  the  difference  between  the 
worst  and  the  moving  average  (which  does  not  use  program  <*  ta) . 
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ABSOLUTE  FDRF.rAST  ERROR  OF  NORMALIZED  DEMAKD 


The  regression-of-slope  algorithm  REGORG)  had  the  smallest  MAD  and  thus 
was  selected  to  be  Included  in  the  linn:  testing  phase.  The  partial-effect 
ulgor Ithin  witli  various  values  of  Che  paiamenter  did  not  do  consistently 
well.  However,  it  d'd  perform  well  in  o study  on  rifle  parts  [GAJDALO  1969 ]. 
It  was  decided  to  Include  this  algorithm  (PARTEF)  with  a parameter  of  .5  so  as 
to  be  as  different  ns  possible  from  the  extremes  of  SEV94  and  MOVAVG  also 
included.  Thus  four  algorithms  were  submitted  for  final  testing  with  the 
Simulator. 
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5.2  Simulation  Results 


In  section  4.2  It  is  stated  that  the  results  of  simulating  each 
algorithm  are  displayed  as  a curve  on  a graph  of  average  requistlon  delay 
(due  to  stockouts)  versus  total  variable  operating  cost.  The  curve  for  the 
best  algorithm  will  be  to  the  left  of  these  for  all  the  others*  indicating  the 
lowest  cost  for  a given  delay  time. 

The  first  graph  shows  the  results  for  the  simulation  run  of  all  items. 

On  this  graph  it  can  be  seen  that  the  SLV94  algorithm  is  much  better  than  a 
simple  moving  average  (MOVAVC) . The  PAI'TEF  algorithm  performed  about  as 
well  as  SEV94.  REGORG,  however,  average!  $130  per  year  per  item  cheaper 
than  SEV94  for  the  current  17-day  mean  rcqulstion  delay.  This  is  a saving 
of  over  52  or  more  than  one  million  dollars  per  year.  It  can  be  seen 
that  the  ranking  of  the  algorithms  studied  In  the  Simulation  is  essentially 
the  same  as  it  was  in  the  screening  phase. 

When  we  consider  Just  the  dynamic  items  in  the  next  graph,  we  see 
that  the  results  are  similar  to  those  for  all  items  together,  except  that 
now  the  differences  between  algorithms  are  more  pronounced.  For  this 
group  of  items  PARTEF  is  noticeably  better  than  5EV94,  but  REGORG  is  still 
the  beet. 
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The  la»t  graph  shows  the  results  l<  • the  rest  of  the  items.  Again  the 
algorithms  that  utilize  program  data  (kl'OKG,  PaRTEF,  and  SEV94)  appear  to  be 
superior  to  the  one  that  doesn't  (MOVAVt  >.  However,  it  cannot  be  definitely 
proven  from  these  results  whether  progran  factors  are  really  more  effective 
than  simple  demand  history  forecasting  algorithms  for  nondynamic  Items  because 
of  the  contamination  with  items  that  art  possible  dynamic.  What  can  be 
determined  is  that  there  is  much  less  difference  among  the  various  program 
factor  algorithms  for  this  group  of  itema  than  for  the  dynamic  items. 

Caution  must  be  exercised  in  using  these  graphs.  While  there  ia  con- 
siderable safety  in  ujlng  them  to  rank  the  algorithms,  there  is  much  less 
assurance  that  the  actual  cost  values  are  correct.  Evaluation  of  these 
costs  is  beyond  the  scope  of  the  present  work. 
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5.3  Lag.  Base  Period  and  Stratification 


The  i Ir.al  tiuce  )>tnphs  Rive  tlx-  r«  ults  of  the  study  to  determine 
the  effects  of  tliet<e  modifications  to  the  forecast  algorithms:  lag 

(Model  4).  Increased  base  period,  and  a stratification  of  the  items  into 
low  and  high  demand  (not  dollar  value)  classes.  This  work  was  done 
with  the  technique  used  for  screening;  thus  rankings  of  many  algorithms 
are  compared.  On  each  graph  is  plotted  the  mean  absolute  forecast  error 
with  the  proposed  modification  versus  that  without  it. 

An  examination  of  the  three  graphs  suggests  the  following  conclusions. 
A two-quarter  lag  degrades  forecasting  by  most  algorithms,  as  does  a 
3-year  base  period.  The  demand-rate  stratification  has  no  overall  effect. 
Stratifications  were  also  attempted  on  unit  price  and  demand  dollar  value 
with  no  effect.  These  are  not  shown  here. 
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The  negative  result  for  the  longer  oaso  period  can  he  challenged 
because  the  advantage  of  the  longer  bas<  period  (smoothing  out  the 
demand)  was  anticipated  by  the  aggregation  into  weapon  systems,  '.'nrk  by 
the  Air  Fence  Logistics  Command  suggests  that  longer  l a so.  periods  are  r.orc* 
effective  with  a moving-average  forecast.  While  this  does  not  seen  likely 
with  the  Army’s  data,  because  ol  the  rr< nds  present,  a longer  bate  period 
may  well  improve  the  forecasts  developc*  with  a program  factor.  .VTfKO 
is  examining  this  possibility  in  a different  project. 

5.4  Final  Conclusion 

The  Army’s  present  program  factor  is  better  than  a simple  nov.ng 
average. 

Regression  on  an  eight-quarter  bast  period  to  determine  the  relation- 
ship between  demand  and  program  (assuming  proportionality)  is  bettor  than  the 
present  foiecasting  method. 

The  savings  on  a ignited  group  of  dynamic  : terns  using  the  recession 
algorithm  are  so  large  tlrut  this  aigcri'htu  car.  be  used  ir  forecast  ing  all 
lteme,  thus  avoiding  the  use  ol  separate  algorithms  e or  aifferent  groups 
of  Items.  In  fact  isolation  of  dynamic  items  from  the  ton -dynamic  was  not 
complete  enough  to  permit,  determination  of  a forecasting  algorithm  preferred 
for  each  group  separately;  nor  does  it  : eem  to  be  cost  effective  to  try  to 
find  separate  algorithms  because  of  the  small  economic  impact  of  the  non- 
dynamic items. 

Lagging  the  demand  as  much  as  two  quarters  does  not  improve  telecasting. 

l.'o  method  was  found  for  classifying  items  based  on  unit  prio.  , average 
demard  rate  or  dollar  \alu«  of  demand  such  chat  different  forecasting  al- 
gorithms would  bo  preferred  for  different  groups  of  items. 
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APPENDIX  1 


CHANGES  TO  SIMULATOR  FOR  PROGRAM  FACTOR  PROJECT 

PREFACE 

The  simulator  used  for  final  selection  of  the  forecasting  algorithm 
is  described  in  the  report,  "ALPHA  4140.39  Simulator"  Hay  1973  by  Martin 
Cohen  (AD  762348).  Changes  to  adapt  It  to  the  present  application,  for 
general  improvements  and  to  correct  errors,  are  described  here  and  keyed 
to  section  headings  in  the  simulator  report  and  to  the  figure  in  the 
Model  Validation  and  Selection  chapter  of  this  report. 

Most  of  the  changes  vara  made  in  order  to  be  consistent  with  the 
recommendations  for  CCSS  in  the  report,  "Evaluation  of  Several  VSL/EOQ 
Models,"  by  Robert  L.  Deemsr  and  U.  Karl  Kruse,  May  1974,  AD  781948. 


INVENTORY  DECISION  PARAMETERS  (SECTION  . .4)  - POD I 41*0.39  INVENTORY 

MANAGEMENT  MOi'EL 


REORDER  POINT 

Instead  ol  using  the  model  111  i.afcty  level  as  defineu  in  "More 
Ado  About  EO(j"  by  V.  J.  i'resutti  and  R.  Trapp  In  NRLQ  17,  12,  June  1970, 
ve  are  now  using  model  IV.  Ip  this  formulation  the  holding  cost  is 
proportional  to  the  sum  of  on-hand  plus  on-order  stock.  The  only  change 
is  to  delete  the  holding  cost  rate,  Cj  , from  the  denominator  of  formula 
(2)  on  page  19  of  the  Simulator  report. 

Another  change  relates  to  the  distribution  assumed  for  the  underlying 
demand  process.  When  the  lead  time  demand  forecast  is  greater  than  20  units 
Presuttl's  Laplace  is  used  as  before,  however,  when  it  is  less  than 
or  equal  to  20,  the  negative  binomial  is  used.  The  theory  is  explained  In 
Deemcr  and  Kruse's  May  1974  report.  Ar.  approximation  is  used  for  the 
negative  binomial;  it  l:  based  on  the  Cnmp-I*aulson  approximation  described 
in  a paper  by  2.J.  n.-rtko,  "Approximating  the  Negative  Binomial,"  in 
Technometrics  b,2,  1966.  See  "Application  of  Negative  Binomial  to  Inventory 
Control,"  ty  K.L.  Deentr,  A.J.  Kaplan  and  W.K.  Kruse,  December  1974, 

AD  A003225,  pages  9-12  tui  the  details. 

The  variance  estinaticn  formula  has  also  been  changed  to  agree  with 
recent  work.  It  is  now  based  on  an  empirical  study  with  AVSCOM  demand 
data,  "Estimation  of  Demand  Variability  Parameters"  by  Alan  J.  Kaplan, 

May  1974,  AL>  731942.  The  final  tables  end  formulas  are  given  in  Deemer 
and  Kruse's  May  1974  report. 

REORDER  QUANTITY 

Formula  (1)  for  0 on  page  14  is  still  correct;  however,  there  are  two 
changes  in  the  way  the  i ir.al  quantity  oi  Jored  is  found. 

Firm,  the  boundary  value  is  no  longer  considered  as  an  alternative 
to  the  resu.1  C ol  the  formula  for  0.  Th.s  means  that  the  so-called  boundary 
opt  Imitation  has  been  abandoned.  The  s«.  *ond  change  is  that  the  actual 
quantity  ordered  includes  the  reorder  point  deficit , the  amount  by  which 
the  assets  were  below  t lie  reorder  point  when  the  buy  was  triggered.  Thus 
the  policy  is  now  to  order  up  to  the  su’  of  the  reotdet  point  and  order 

4o 


quantity,  which  Is  more  in  line  with  actual  practice  and  with  the  assump- 
tions on  which  the  inventory  model  is  based. 

Boundary  optimization  has  been  abandoned  because  now  the  actual 
quantity  ordered  Includes  the  deficit,  which  would  bring  it  above  the 
boundary  at  which  it  was  set  to  minimize  cost.  Rather  than  adopt  a more 
complicated  procedure  to  try  to  save  a questionable  policy  like  boundary 
optimization,  the  whole  Idea  has  been  dropped  under  CCSS. 

To  summarize  the  changes  in  the  reorder  quantity,  the  former  policy 
(pages  17  and  18  of  the  May  1973  report)  was  to  find  the  smallest  feasible 
reorder  quantity  that  minimized  the  projected  Inventory  related  variable 
cost,  and  to  order  exactly  that  quantity  (except  when  the  deficit  was  too 
large,  in  which  case  a multiple  of  that  quantity  was  ordered).  The  new 
policy  is  to  order  the  deficit  as  well  as  the  order  quantity  and  not  attempt 
to  minimize  the  projected  cost  (except  the  minimization  Implied  in  the  Q 
formula  itself)  by  considering  the  Q value  at  the  order  cost  boundary. 

FREQUENCY  OF  RECOMPUTATION  OF  REORDER  POINT  AND  QUANTITY 

On  page  13  of  the  Simulator  report  it  is  stated  that  the  reorder  point 
and  reorder  quantity  are  recomputed  each  time  the  reorder  point  is  creased 
as  well  as  once  a yaar  for  LDV  items  and  monthly  for  MDV/HDV  items.  This 
has  been  changed  with  the  installation  of  a "buy-back"  policy.  Nov  computa 
tions  are  not  made  when  the  reorder  point  is  crossed,  and  the  recomputation 
frequencies  have  been  changed  to  semi-annually  for  LDV  items  and  quarterly 
for  MDV/HDV  items.  Also,  reparable  items  are  now  considered  to  be 
MDV/HDV  items,  and  the  item's  class  (LDV  or  MDV/HDV)  is  now  determined 
each  time  the  decision  is  to  be  made , using  the  current  value  of  the  demand 
forecast.  The  buy-back  policy  assumes  a 15  day  reduction  in  procurement 
lead  time;  this  has  been  implemented  in  the  simulation  except  in  the  case 
that  a buy  actually  occurs  at  the  same  time  as  a recomputation.  For  more 
information  see  "Frequency  of  Requirements  Determination"  by  Robert  L. 
Deemer,  October  1974,  AD  A003227. 
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ZERO  DEMAND  FORECAST 


On  page  18  of  the  Simulator  report  t is  stated  that  the  safety  level 
component  of  the  reorder  point  is  set  to  tero  and  the  reorder  quantity  is 
set  to  one  if  Che  demand  *ate  or  varietur  forecast  is  zero.  This  has 
been  changed  because  it  to^y  be  very  unrealistic  for  LDV  items.  Now,  the 
number  of  months  in  the  safety  level  and  reorder  quantity  are  calculated 
with  forecasts  that  are  not  permitted  to  be  less  than  one-third  unit. 

The  number  of  units  in  these  levels  are  then  calculated  by  multiplying  the 
number  of  months  times  the  demand  rate  forecast  (which  may  be  less  than 
one-third  unit  per  year  or  one-thirty-si::th  unit  per  month).  Whenever 
this  is  done  the  resulting  reorder  quantity  is  forced  to  be  at  least  one 
unit;  this  is  done  to  avoid  placing  several  orders  for  a fractional  number 
of  units. 
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IUI1 IALIZAT10M  ANiJ  'VARML'P  (SECTION  2.5) 

Before  simulating  each  item  an  lit j r a]  demand  forecast  is  made  and 
start-up  value-:  are  computed.  Then  the  item  is  simulated  ter  some  time 
without  collecting  performance  statistics  in  order  to  reduce  the  effects 
of  the  assumptions  made  in  estimating  the  initial  values.  It  has  been 
possible  to  Increase  this  warmup  time  from  one  year  to  two  years  because 
we  now  have  seven  years  of  demand  data.  The  initial  forecast  is  now  based 
on  the  first  two  years  of  demand  rather  than  the  final  two  years  because 
It  was  noticed  that  there  was  a general  down  trend  in  demand  during  the 
later  years  of  the  data.  The  effect  of  having  perfect  information  in 
the  first  two  years  is  risked  by  the  two  year  warmup. 

Since  the  purpose  of  this  project  is  to  compare  different  forecast 
ulgoritluns,  the  algorithms  should  all  have  the  same  starting  conditions, 
i.e.,  the  conditions  when  statistics  are  first  accumulated,  just  after 
warmup.  In  order  to  accomplish  this,  wa-mup  is  now  exactly  the  same  for 
all  the  algorithms.  This  is  done  by  using  the  moving  average  forecast 
algorithm  during  warmup,  no  matter  which  algorithm  is  to  be  used  afterward. 
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OUTPUT  (SECTION  3.3)  ASSET  MONITOR  PERFORMANCE  CALCULATIONS 


REQUISITIONS  SHORT 

Time-weighted  requisitions  short,  from  which  the  days  wait  measure 
is  calculated,  is  now  computed  with  a more  accurate  approximation.  The 
basic  formula  for  time-weighted  units  short  is  still  the  same  as  given  in 
Section  2.C  of  the  Simulator  report.  Now,  however,  a better  method  is 
used  to  estimate  the  number  of  requisitions  represented  by  these  units. 
Rather  than  3imply  estittfating  the  requisitions  short  from  the  units 
short  by  dividing  by  the  current  order  size  as  was  done  in  the  previous 
version,  now  the  number  of  requisitions  short  is  carried  forward  as 
a total  and  updated  by  reducing  or  increasing  its  value  by  the  same  fraction 
as  the  change  in  units  short. 

For  example  suppose  that  at  one  point  in  simulated  time  there  is 
estimated  to  be  2 requisitions  for  a total  of  10  units  of  stock  on 
backorder.  Assume  that  two  more  requisitions  are  placed  for  a total  of 
six  units  and  that  five  units  are  received  from  the  supplier.  After 
this  there  will  bel0+6-5*tll  units  backordered;  hut  how  many 
requisitions?  The  old  way  would  have  est imated  the  number  of  requisitions 
by  dividing  the  number  of  units  by  the  order  size  in  the  current  period, 
i.e.  6/2  * 3 units,  then  requisitions  short  ■ 11/3  * 3.A7.  The  new  method 
estimates  the  additional  requisitions  on  backorder  as  the  same  percentage 
increase  as  the  units  (11-10) /10  « 10%;  then  requisitions  short  * 

2 + 10%  * 2.2  in  the  example. 

The  new  values  have  tended  to  be  about  one-third  less  than  the  old 
one  and  are  considered  much  more  reasonable. 

OBSOLESCENCE  FENALTY 

Instead  of  including  ar.  obsolescence  rate  in  with  the  holding  cost 
rate,  now  an  obsolescence  penalty  is  computed  directly  and  added  to  the 
holding  cost.  This  penalty  is  computed  as  follows.  At  the.  end  of  the 
simulation  of  each  item,  any  on-hand  or  cn-order  stock  (assets)  in  excess 
of  the  sum  of  the  reorder  point  ar.d  reorder  quantity  is  reduced  by  the 
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present  value  of  total  future  demand  and  multiplied  by  the  item's  unit 
price.  To  project  future  demand,  the  ptogram  data  for  the  period 
following  the  simulation  is  set  to  the  .iverage  of  the  values  for  the  last 
two  years.  The  rationale  for  this  method  of  estimating  obsolescence  is 
explained  in  the  Deemer  and  Kruse  May  1974  report  on  pages  15-17. 
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APPENDIX  2 


AGGREGATE  NORMALIZED  DEMAND  QUANTITY  AND 
FLYING  HOURS  BY  AIRCRAFT  SYSTEM 
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